Calculation of magnetic exchange couplings in S=3/2 honeycomb system 
Bi3Mn40i2(N03) from first principles 
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Absence of magnetic ordering in Bi3Mn40i2(N03), (BMNO) which lias a magnetic subsystem 
that consists of honeycomb bi-layers of Mn^^ ions with spin S=3/2, has raised the expectation that 
its ground state is strongly frustrated due to longer-range antiferromagnetic interactions. We cal- 
culate the magnetic exchange coupling constants of the BMNO complex within a density functional 
approach and find that the dominating interactions between Mn spins are the antiferromagnetic 
nearest-neighbor Ji and interlayer interaction Jc. The largest interaction is Jc, which substantially 
exceeds Ji. Longer-range interactions are antiferromagnetic, but not strongly frustrating. 
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In frustrated magnets spins interact through compet- 
ing exchange interactions that cannot be simultaneously 
satisfied, giving rise to a large degeneracy of the sys- 
tem's ground state. The resulting fluctuations of the 
spins are at the root of remarkable collective phenomena 
such as emergent gauge fields and fractional particle exci- 
tationsi. In this context the particular magnetic proper- 
ties of Bi3Mn40i2(N03) (BMNO) have recently received 
considerable attention^r^i. BMNO is a layered manganese 
oxide with as its main magnetic building block a honey- 
comb lattice of Mn^"*" ions, each carrying a spin S=3/2. 
These manganese honeycomb lattices form bilayers (see 
Fig. 1), with interactions between the spins that are an- 
tiferromagnetic. In magnetic susceptibility and specific 
heat studies the absence of long-range ordering down to 
0.4 K was reported'^, which was confirmed in a high-field 
ESR study-^. It was subsequently observed by inelastic 
neutron scattering^ that long-range spin order can be in- 
duced by an external magnetic field. 

The lack of spontaneous magnetic ordering within the 
Mn honeycomb lattices down to very low temperature 
was proposed to be due to frustration of magnetic in- 
teractions, in particular a competition between first and 
second neighbor antiferromagnetic Mn-Mn exchange in- 
teractions^. Currently BMNO is therefore viewed as a 
model system with strongly fluctuating S=3/2 spins on 
honeycomb lattice, with a magnetic frustration that is 
due to the presence of longer-range antiferromagnetic ex- 
change interaction&2r— . 

The usefulness of this paradigm critically depends on 
the exact values of the exchange magnetic interactions 
in this material. The exchange coupling constants have 
been extracted from inelastic neutron scattering experi- 
ments^, where however the values of the couplings are not 
fully constrained. Theoretically, coupling constants have 
been estimated from unrestricted Hartree-Fock calcula- 
tions and a perturbation method^, which puts forward 
that a ferromagnetic third nearest neighbor interaction 
is responsible for the absence of long range magnetic or- 
der in BMNO. We performed ab initio calculations of the 
magnetic exchange coupling constants of BMNO based 
on density functional theory and find that the nearest- 



neighbor Ji and interlayer interaction Jc are the most 
important interactions. In-plane longer-range interac- 
tions, including the third neighbor one, are antiferromag- 
netic and much smaller. Taking all these into account we 
find the system is significantly away from full frustration 
within the honeycomb planes. In contrast to what has 
been asserted so far for BMNO, we find that the domi- 
nant interaction is the coupling between the honeycomb 
bilayers, Jc, with a value that can be as large as twice Ji. 
The existence of interlayer interaction was recognized in 
a recent study by Ganesh et al. in Ref . , but our find- 
ing that it actually dominates the exchange requires an 
alternate perspective on the factual origin of the observed 
magnetic frustration in BMNO. 

We extract the exchange coupling constants in BMNO 
within GGA-I-U for different values of the Coulomb repul- 
sion U on the Mn site. So far BMNO has not been stud- 
ied with first-principles calculations and we use frozen 
coUinear spin conflgurations to calculate very accurately 
the energies for different spin configurations. For the 
calculations we used an all-electron full potential local- 
orbital scheme implemented in FPLO code^'^° with basis 
sets: Bi (5s, 5p, 5d, 6s, 6d, 6p, 7s, 7p), Mn (3s, 3p, 3d, 
4s, 4p, 4d, 5s), N (Is, 2s, 2p, 3s, 3p, 3d), and O (Is, 2s, 
2p, 3s, 3p, 3d) states. The exchange-correlation energy 
functional was evaluated within the generalized gradi- 
ent approximation (GGA) using the Perdew, Burke, and 
Ernzerhof parameterizationii. The number of irreducible 
k points used in the calculations were 48 from a sampling 
of 4x4x3 mesh. 

We use experimental crystal structure of 
Bi3Mn40i2(N03), which is trigonal with space group 
P3 and consists of two MnOg layers separated by Bi 
and NO3 layers (see figure [Tp. The edge-shared MnOe 
octahedra lie in ab plane and form honeycomb lattice as 
shown in figure [Ijb) . The four inequivalent Mn atoms 
form two hexagonal planes separated by Bi atoms. 
These hexagonal planes are in distorted environments of 
MnOg octahedra resulting in small differences between 
the planes. We simplify the situation by treating 
the two layers as if they were the same, thus having 
the same in-plane magnetic interactions. The main 
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FIG. 1: (color online). Crystal structure of BMNO. (a). Dis- 
torted MnOe octahedra are edge sharing with oxygen at the 
vertices, forming hexagonal bi-layers. For clarity the bottom 
hexagonal layer is repeated along c direction, which empha- 
sizes the hi- layer structure, (b). Top view of the network 
of edge-sharing octahedra in BMNO. The honeycomb lattice 
is indicated by the dark black lines that connect Mn atoms 
(black dots). The oxygen atoms shared by the MnOs units 
form layers. The orange (medium gray) spheres are Bi atoms, 
the yellow (light gray) spheres nitrogen atoms, and red (dark 
gray) spheres are oxygen atoms. 



purpose of presenting this streamlined picture is to 
understand and illustrate nature and main trends of 
magnetic interactions. But the small crystallographic 
differences will cause changes in the details of magnetic 
interactions that are presented here. Five different 
magnetic spin configurations (as shown in Fig. [2]) are 
considered to calculate the magnetic interactions within 
the honeycomb lattice and between the bilayers with 
the GGA+U method for U = 4.0, 6.0, and 8.0 eV at the 
transition metal Mn site, which is within the range that 
of U values extracted from spectroscopic measurements 
(4-5 eV)—^ and used in other ah initio calculations 
(4-10 eV)ii-il. The Hund's rule couphng Jh is 1 eV. 

Total energies for different Mn spin configurations in 
a supercell with 92 atoms were obtained in the GGA+U 
formalism employing the atomic limit functional. The 
calculated relative energies for the spin configurations A, 
B, C, D, and E with the different U values are listed in 
Table [H where the energy of the most stable configura- 
tion "D" energy was used as the zero of energy. With 
the energies of five ordered spin states, all four exchange 



FIG. 2: (color online). Schematic representation of five or- 
dered spin states A, B, C, D, and E in a 2x2x1 supercell 
of BMNO used for the determination of exchange constants. 
Only Mn atoms with the spin exchange paths are shown. Red 
filled (dark gray) circles are used to highlight spins that are 
flipped with respect to the reference configuration. Three in- 
plane magnetic interactions given as Ji, J2, and J3, whereas 
Jc represents the inter-bilayer interaction. 





(7=4.0 eV 


17=6.0 eV 


(7=8.0 eV 


A 


0.3745 


0.2213 


0.1118 


B 


0.1180 


0.0684 


0.0335 


C 


0.1687 


0.1013 


0.0518 


D 











E 


0.2983 


0.1752 


0.0852 



TABLE I: Relative energies in eV of different spin-states 
structures (as shown in Fig. [2]) of BMNO obtained from 
GGA+U {U = 4.0, 6.0, and 8.0 eV) calculations. 



constants ( Ji, J2, and J3, between first, second and third 
nearest neighbors, respectively and Jc the inter-bilayer 
coupling) can be determined. 

To evaluate the magnetic interactions, we map the 
DFT energy differences of five magnetic spin configura- 
tions onto a classical spin-3/2 Heisenberg model. The 
magnetic exchange coupling constants are calculated 
from the following equations: 
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In Table In] Ji, J2, J3, and Jc are listed as calculated 
from the relative energies of these spin configurations, 
for different values of U . All interactions come out to be 
antiferromagnetic and we observe that the dependence 
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exchange 


U = 4.0 eV 


U = 6.0 eV 


U = 8.0 eV 


Ji 


3.00 (34.8) 


1.7 (19.7) 


0.9 (10.4) 


J2 


0.4 (4.6) 


0.2 (2.3) 


0.1 (1.2) 


Ja 


0.47 (5.4) 


0.3 (2.7) 


0.2 (2.3) 


Jc 


4.10 (47.6) 


2.7 (31.3) 


2.1 (24.4) 


J2/J1 


0.13333 


0.1176 


0.1111 


■h/Ji 


0.15667 


0.17647 


0.2222 


Jc (direct) 


2.97 


1.88 


1.08 



TABLE II: Exchange constants in BMNO in meV calculated 
by FPLO GGA+U with the atomic limit (AL) functional, on 
the basis of Eqs.[T]|4l Slater parameters are chosen as Fo = U , 
F2 =8.6 eV and F4 = 5.4 eV, i.e. J = (F2 + F4)/14 = 1 eV. 
Values in Kelvin (K) in parentheses. The direct evaluation of 
Jc is based on the energy difference between ferro and anti- 
ferro coupled FM honeycomb planes. 

on U is rather moderate as far as the relative values of 
the exchange constants is concerned: U roughly tends 
to scale the different exchange interactions in the same 
manner. The in-plane Ji is almost an order of magnitude 
larger than J2. It is also considerably larger than J3, so 
that these interactions are only weakly frustrating. The 
magnitude of the couplings J2 and J3 are similar. As 
the bond length corresponding to Ji (2.8697 A) is much 
shorter than that for J2 or J3 (4.9692 A/ or 5.7383 A), 
it is not surprising that Ji is the strongest interaction 
vi^ithin the honeycomb plane. However, the largest ex- 
change interaction is the interlayer Jc, which is also an- 
tiferromagnetic and has about two times the value of Ji 
for larger values of U. 

Let us compare these results with the exchange con- 
stants found by inelastic neutron scattering. Matsuda et 
al. have investigated magnetic ordered and disordered 
ground state using inelastic neutron scattering experi- 
ments^. From the data they extract a constraint that 
relates the spin-spin interactions based to the observed 
magnetic excitations: 3 Ji + 6 J2 -I- Jc ~ 6 meV. By subse- 
quently using in addition an approximated value of Jc/ Ji 
= 1/2 from their experiments and the critical value of 
J2/ Ji = 0.15 from literature^^^, which ensure very strong 
frustration and are consistent with the observed absence 
of long-range magnetic ordering, Matsuda et al. obtain 
Ji = 2 Jc = 1.4 meV and J2 = 0.20 meV. Compared to the 
experimental analysis, we consider in our theoretical cal- 
culations a fourth interaction (J3) corresponding to the 
third nearest neighbors interaction within the honeycomb 
lattice. In this case the experimental constraint of Mat- 
suda et al. is C=3 Ji -I- 6 J2 + Jc + 3 J3 ~ 6 meV, where 
we find from our calculations C=5.9 meV, for U = 8.0 eV. 
The calculated constraint and therefore also the calcu- 
lated exchange interaction are in full agreement with the 
dispersion of magnon excitations that were measured by 
inelastic neutron scattering. We obtain a ratio of J1/J2 



that is similar to the one that was used by Matsuda et 
al.^ . However, our interlayer coupling Jc is larger than 
the in-plane coupling constants, contrary what was calcu- 
lated from inelastic neutron scattering experiments^. To 
cross-check this result, we also calculated the energies of 
fully ferromagnetically polarized honeycomb planes that 
are coupled either ferro- or antiferromagnetically. This 
allows for the direct evaluation of the intra-plane ex- 
change but has the disadvantage that because a state 
with ferromagetic planes is the one furthest away from 
the in-plane antiferromagnetic ground state, significant 
renomalizations of the inter-plane magnetic bonds can 
occur. As shown in Table HIl we find in this case a ratio 
of Jc/ Ji closer to one, the exact value of which depends 
on U. This is a substantial renormalization, but does 
not change the conclusion that the inter-layer coupling is 
the leading exchange interaction. This relatively large Jc 
is related to the fact that the in-plane Mn-O-Mn angles 
are close to 90'', implying a strong competition between 
antiferromagnetic and ferromagnetic superexchange, of 
which it is a priori not clear which will win. The inter- 
plane exchange paths, via the covalent bismuth states, 
do not suffer from such a geometric frustration of inter- 
actions. 

Thus from the theoretical evaluation of magnetic ex- 
change coupling constants for BMNO, we find all of them 
(first-neighbors, second, and as well as third neighbors 
within the honeycomb plane and between the planes) to 
be antiferromagnetic. The calculated values fulfill the 
constraint that is put on the weighted sum of exchange 
interactions by inelastic neutron scattering datct^. In 
this constellation, only the second nearest neighbor in- 
teraction J2 would be frustrating a full-blown antifer- 
romagnetic ordering. As J2 is comparable to J3, both 
of which are much smaller than Ji, which is again ex- 
ceeded by Jc, we find from the calculations that the ef- 
fect of this frustration is not strong. We conclude that 
the observed absence of magnetic ordering down to tem- 
peratures lower than the temperature scale set by the 
calculated exchange interactions requires an explanation 
that goes beyond in-plane magnetic frustration. Possi- 
bly oxygen vacancies are of importance which will not 
only lead to disorder of magnetic bonds but might also 
lead to doping of the manganese layers, in which case 
the presence of Mn^+ would lead to a competition be- 
tween ferromagnetic double exchange and antiferomag- 
netic superexchange, which is known to strongly affect 
the magnetic properties and ordering in weakly doped 
Mn'^^ perovskite manganites-?. At any rate, a descrip- 
tion of the magnetic properties of BMNO will have to 
considered the dominating coupling constant, inter bi- 
layer exchange, explicitly. 

We thank A. Paramekanti for stimulating our interest 
and discussing this honeycomb manganite with us and 
A. Tsirlin for discussions. 
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